Re-entrant spin susceptibility of a superconducting grain 
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Abstract 



We study the spin susceptibility x of a small, isolated superconducting grain. 
Due to the interplay between parity effects and pairing correlations, the de- 
pendence of x on temperature T is qualitatively different from the standard 
BCS result valid in the bulk limit. If the number of electrons on the grain 
is odd, x shows a re-entrant behavior as a function of temperature. This 
behavior persists even in the case of ultrasmall grains where the mean level 
spacing is much larger than the BCS gap. If the number of electrons is even, 
x{T) is exponentially small at low temperatures. 
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By now it is well-known that the properties of an isolated, mesoscopic superconducting 
grain are quite different from those of a bulk sample!. First of all, since such a grain 
carries a fixed number, N, of electrons, its behavior depends strongly on whether N is even 
or odd. Second, fluctuation effects become important as the size of the grain decreases. 
The interplay between parity and fluctuation effects crucially depends on the ratio S/Aq 
of two characteristic energies: the mean level spacing 5 and the bulk superconducting gap 
Aq. As long as the grain is not too small, 5 <C Aq, the fluctuation region AT around the 



critical temperature T c is narrow, AT/T C ~ yo/An <C 1, and the mean field description 
of superconductivity is appropriate. Parity effectsai appear at temperatures lower than 
a crossover temperature T e g- ~ A / In J 87tAq/ 5 2 which, in the experiments!, is typically 
of the order of 10 — 30% of T c . The dependence of T eS on A signals that the even-odd 
asymmetry is a collective effect due to pairing correlations. As the size of the grain is 
decreased, fluctuations start to smear the superconducting transition!. The finite level 
spacing suppresses the BCS gap in a parity-dependent wayoQ. When 5 becomes of the order 
of A , AT ~ T c and the BCS description of superconductivity breaks down even at zero 
temperature^. The regime 5 > A is dominated by strong pairing fluctuationsBTLj. 

The present-day interest in ultrasmall superconducting grains was triggered by the exper- 
iments of Ralph, Black and TinkhamB, who were able to contact a single, nanometer-sized 
Al grain with current and voltage probes. They obtained tunneling spectra that revealed 
the presence of a parity-dependent spectroscopic gap, larger than the average level spacing, 
which could be driven to zero by an applied magnetic field. In this Letter we propose to 
measure the temperature dependence of a thermodynamic quantity - the spin susceptibility 
X - as a means to detect both parity effects and pairing correlations. As we will show be- 
low, pairing correlations give rise to a specific temperature dependence of thermodynamic 
quantities. This enables a more quantitative investigation of fluctuation effects^!. 

Spin paramagnetism of small particles has been considered in the pastel and very recently 
parity effects in the susceptibility were measured for an ensemble of small, normal metallic 
grainsE-3. Spin susceptibility is very sensitive to BCS pairing as well. Yosidall! showed that, 
due to the opening of the superconducting gap, x vanishes at zero temperature^. We will 
show that the combined effect of parity and pairing introduces qualitatively new features in 
the temperature dependence of x- Most interestingly, these effects might be observed even 
in the regime 5 ^ Ao- The results of our work are summarized in Figs. |T] - |3|, where we 
plot x as a function of temperature T for odd and even parity. In particular we want to 
emphasize that the odd susceptibility shows a re- entrant behavior as a function of T for any 
value of the ratio S/A . This re-entrance is absent in normal metallic grains; it is a genuine 
feature of the interplay between pairing correlations and parity effects. 

The BCS pairing Hamiltonian for a small grain can be written as!H! 

W= ( 6 « _ WbH) C ] n u C n ^ - X5 B L B n, (i) 

n,cr=± m,n 

where = cl n>+ cl n _. The indices m, n label the single particle energy levels with energy e m 
and annihilation operator c m (7 . The quantum number a — ± labels time-reversed equally 
spaced states (with an average spacing ~ 5 = 1/Vo, where vq is the density of states at 
the Fermi energy). The external magnetic field H couples to the electrons via the Zeeman 
term, \i& is the Bohr magneton. We put the g-factor equal to two, ignoring any spin- 
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orbit effects (see Ref. ^0|). At the low magnetic fields of interest here, we can neglect the 
orbital contribution to magnetic energy, as it is smaller than the Zeeman energy by a factor 
~ (kpr)(Hr 2 /$o) ( r is the size of the grain and $o the flux quantum). Finally A is the 
dimensionless BCS coupling constant. Since the Hamiltonian contains only pairing terms, 
an electron in a singly occupied level cannot interact with the other electrons. 

The spin susceptibility of a grain with an even (e) or an odd (o) number N of electrons 
is defined as 



Xe/o(T) = - 



d 2 f e/o (T,H) 



OH 2 



(2) 



H=0 



where jF e / Q = — ThiZ e / is the free energy of the grain and the partition function Z(T,N) 
should be evaluated in the canonical ensemble. We will perform the calculation with the 
help of a parity projection technique!!! and by means of exact canonical methods based 
on Richardson's solution!. The grand partition function reads 



Z e/o (T,fi) 



(1/2) £ e^ T [l±e™ N 

N=0 



}Z(T,N) = (1/2) (Z + ±Z_). 



(3) 



The partition function Z + is the usual grand partition function at temperature T and chem- 
ical potential /i + = fi. The grand partition function Z_ describes an auxiliary ensemble at 
temperature T and chemical potential /i_ = [i + iirT; it is a formal tool, necessary to include 
parity effects. The chemical potential /i will be placed between the topmost occupied level 
and the lowest unoccupied level in the even case, while it will be at the singly occupied 
level in the odd case. Since we are interested in the evaluation of fluctuations effects, it is 
convenient to express the grand partition functions Z± using the path integral formulation 
of superconductivityBIlll, 



(4) 



/£> 2 Aexp< 




Trl n (l-G»A)-^; 


} 


/ V 2 A exp < 




f 



Here, /3 = 1/T and Z\ is the partition function for non-interacting electrons. The matrix 

Green function G\ is given by G^{e u ) = {iuj v + fi B H)cr^ — (e n — fM±)(T^ , where to v is 

a fermionic Matsubara frequency, [i = x, y, z) are the Pauli matrices, and is the 
identity. Finally, the matrix A is given by A = (A/2)(a x + ia y ) + h.c. A direct calculation 
of the partition function (f|) is impossible in general. Below, we first discuss two limiting 
cases which are tractable analytically: S/Aq <C 1 and S/Aq ^> 1. Then we present the 
complete temperature dependence of the spin susceptibility evaluating Eq. (ffl) numerically 
for arbitrary values of S/A , with the help of the static path approximationa 

Large grains ( A 3> 5) . In this limit it is sufficient to evaluate the partition function 
in a saddle point approximation, since fluctuations will not contribute significant ly@. As a 
result we find 



Xe/o 



cosh 



-2 £ 4 



2T 



=F Z_ sinh 



-2 E. 



2T 



z,±z_ 



(5) 
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where E± >n = ye^ + A±. The saddle point values of A± are found from the equations 

A=?4Ei; th I 2T J' (6) 



n,<r 



The partition functions for the two ensembles are 



exp<^ 



In 2 




AT 
\5T 



(7) 



where E± n = E± >n — cr^sH and £ n = e n — ji. 

At low temperatures T A , the ratio Z_/Z + can be calculated easily; one finds 



Z-/Z + ~ 1 — 8itT A 1 5 2 exp (—(3A ). Parity effects are important if this ratio is ~ 1, i.e. 
at temperatures T < T cS . At temperatures T eff C T C Ao, parity effects can be ignored 
and the spin susceptibility is found to decrease exponentially, as in the BCS case, 



For T <C T e g, Eq. ([5]) can be approximated as 

_ 87T/4Ao 2gAn _ /4 

Ae — ^2 ' ^° — y ' 

We see that Xe remains exponentially small, like in the BCS case (||), but with an exponent 
— 2/3 Ao rather than — (3Aq. This reflects the fact that excitations are actually created in 
pairs. In odd grains, the unpaired spin gives rise to an extra paramagnetic (Curie-like) 
contribution to the spin susceptibility. As a result Xo will show a re-entrant effect at low 
temperatures (see Fig. 0). Although the re-entrant behavior is essentially a single electron 
effect, we stress that it can be detected experimentally using granular systems with many 
well-separated grains (to avoid collective effects due to tunneling). Such systems contain 
even grains as well, but their susceptibility is exponentially small at the temperatures of 
interest; thus their contribution to the response of the system will be negligible. 

Ultrasmall grains ( Ao -C S) . A reduction of the grain size leads to a suppression of the 
gap A. For ultrasmall grains with A <C 5, the mean field approximation gives A = 0: the 
grain behaves as a normal metal. The non-interacting, parity-dependent spin susceptibility 
can be found from Eq. (H), see the topmost curves in Figs. |2] and [| Note in particular the 
monotonous dependence of \ Q on T. The temperature scale at which parity effects appear 
is set by the average level spacing. If T ^> 5, parity effects are exponentially small and 
Xe/o(T) — Xp[1 T (2T/5) exp(— 7t 2 T/<5)]. In the opposite limit T <C 5, Xe is exponentially 
small, Xe(T) — (8/i^/T)e~ /3<5 , as we need to excite an electron out of the topmost, doubly 
occupied single particle level to magnetize the grain. For an odd grain, Xo{T) — [J>%/T at 
T ^ 5: the topmost level is occupied by a single electron that gives a Curie-like contribution. 

The saddle point approach entirely ignores the fact that the fluctuation region AT 
around T c grows as the size of the grains is reduced. Due to the presence of fluctua- 
tions, the behavior of small grains will be different in a distinct way from normal metal- 
lic grains. In the limit T 3> 5, both fluctuation and parity effects are small; it there- 
fore suffices to consider the fluctuation correction <5xfluc to xp> evaluating Z + , Eq. (§]), 
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in Gaussian approximation. As a result, we find Sx&uc/xp — ~ 5/2T\ia(T/ A ); hence 
Xe/o(T) — Xp[1 T (2^/(5) exp(— ti 2 T/5) — 5/2Tln(T/A )]. Superconducting correlations sup- 
press the susceptibility; due to its algebraic dependence on T this suppression is stronger 
than the parity correction at temperatures T > <51nln<5/A . In the opposite limit, T ^ 5, 
fluctuations are strong and the Gaussian approximation fails. However, the susceptibility 
can still be obtained analytically by considering a few levels close to the Fermi energy with 
a renormalized pairing interaction A = l/ln(5/Ao)&0. Consider first a grain with an even 
number of electrons. It costs an energy ~ 5 + 5/ln(5/A ) to excite an electron from the 
topmost, doubly occupied level to the lowest unoccupied level. Correspondingly, the leading 
temperature dependence of the spin susceptibility is 

Xe (T) ~ ^e-^i+m-V/Ao)) + Q^-WSy (9 ) 

The even susceptibility is exponentially small, like in the case of a normal metallic grain, 
but with an exponent —/35(l + \), rather than —(35. Similarly, we find the spin susceptibility 
for a grain with an odd number of electrons 

Xo (T) ~ ^[l + 8e -/3a(2+m-V/Ao))]. (10) 

The paramagnetic contribution from the single spin dominates at all temperatures below 5. 
Compared to the case of a normal metallic grain, the odd susceptibility is non-monotonous: 
upon lowering temperature, \ Q first decreases due to superconducting fluctuations; at tem- 
peratures T ~ 5 a re-entrant behavior sets in which persists down to the lowest temperatures. 

Reentrant susceptibility. The various limiting cases discussed so far provide evidence for 
the appearance of an anomaly in the spin susceptibility Xo- For large grains (A ^> 5) the 
mean field approximation, Eq.(|5|), leads to the re-entrant behaviour shown in Fig. |l|. We 
will show that this is a unique signature of pairing correlations which is present even in 
ultrasmall grains. To this end we study the complete temperature dependence of Xe/o for 
arbitrary values of the ratio 5/A . 

The physics of re-entrant susceptibility can be grasped by evaluating Eq. (|j), in the 
static path approximation^. This amounts in retaining only the static fluctuations (beyond 
the Gaussian approximation) in the path integral. In Figs. |2| and |3] we show the results 
of this calculation for the odd and the even cases respectively. The re-entrance in the odd 
case is visible even in systems with a ratio 5/A ~ 50 (!) and provides the signature of the 
existance of pairing correlations in the ultrasmall regime. Frequency dependent fluctuations 
will further enhance the size of the reentrant effect. The results are plotted for a system 
of N = 200 electrons at half filling and the BCS coupling is chosen to fix the ratio 5/ A 
(A ~ 0.1 — 0.2). The merit of the static path approximation combined with the analytic 
analysis in the limiting cases is that it allows to obtain a coherent quantitative physical 
picture in the whole temperature range. 

As a final check of our results we computed Xo using the exact solution of Ref.i. The 
result is presented in Fig.|2| (dashed lines). As expected, the re-entrant effect is slightly larger 
(~ 15%). In order to obtain this result we considered all the different states with excitation 
energy up to a cutoff A ~ 405 for a system with iV < 100 electrons. In the inset we show 
the scaling analysis for different N and different energy cutoffs. This analysis becomes more 



5 



and more difficult upon increasing temperature because of the exponential increase of the 
number of excited states needed. 

In this Letter we proposed the study of the spin susceptibility of a metallic grain as a 
very sensitive probe to detect superconducting correlations. For grains in the nanometer size 
regime the odd spin susceptibility is a unique signature of pairing. In grains of dimensions 



of the order of few nanometers as those studied in Refs. 14 the reentrance should be of the 



order of 10 — 20% of the Pauli value and it could be measured using the techinque used in 



Ref. 17 
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FIGURES 
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FIG. 1. Ao 3> 5 - Spin susceptibility as a function of temperature T/T c for an odd (upper 
curve) and even (lower curve) grain, respectively. The BCS result (middle curve) is reported for 
comparison. The susceptibility is normalized to its bulk high temperature value xp = 2//#/5. We 
used 5/Aq = 0.1. 




12 3 4 

T/8 

FIG. 2. Ao < S - Reentrant spin susceptibility as a function of temperature for an odd 
grain in the static path approximation (solid lines); the top-most curve without re-entrance is 
the non- interacting limit. The dashed lines are obtained by an exact canonical calculation (see 
text), using N = 100 electrons and a large enough maximum excitation energy (A ~ 405). In 
the inset the dependence of the susceptibility as a function of N and A. Curves are labelled with 
(TV, A). 
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FIG. 3. Ao < 5 - Spin susceptibility as a function of temperature for an even grain. As in the 
previous figure, the thin solid line is the non-interacting limit. 
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